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Abstract 

It is reported a combined numerical approach to study the localiza- 
tion properties of the one-dimensional tight-binding model with poten- 
tial modulated along the prime numbers. A localization-delocalization 
transition was found as function of the potential intensity; it is also 
argued that there are delocalized states for any value of the potential 
intensity. 
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Problems related to the distribution of prime numbers are usually inter- 
esting and difficult. Maybe, the most challenging open problem in mathemat- 
ics is the so-called Riemann's hypothesis on the distribution of the (nontriv- 
ial) zeros of his zeta function, a problem directly related to the distribution 
of prime numbers and other mathematical issues. Also, the interplay be- 
tween number theory, with the set of prime numbers playing a major role, 
and physics has been fruitful; see the short note by Wolf and references 
therein for recent discussions. 

In the last years there has been a large interest in spectral properties of 
discrete Schrodinger operators on Z 2 (2) 

(H v u) n = u n+l + M n _! + XV n u n , 

with nonrandom finite valued potentials V = {V n } (see [0, [| ||, 0, [7J and 
references therein); A > is the potential intensity. Such one-dimensional 
models with potentials along almost periodic sequences taking a finite num- 
ber of values 0, [3], |J, like substitution and circle sequences, have been dom- 
inated by singular continuous spectra. Recall, however, that for the cases of 
periodic and random potentials V it is well known that the hamiltonian Hy 
has absolutely continuous and point spectra, respectively ||; these are math- 
ematical characterizations of delocalized and localized (insulator) behaviours, 
respectively. 

There has also been recent work on Schrodinger operators with sparse 



potentials ||, [10], [XT], [l^, |T3|| . A characteristic defining a sparse potential 
is the constancy of its values in a sequence of increasing intervals. In this 
case one can have point, absolutely continuous and singular continuous spec- 
tra, depending on the growing speed of the gaps and how the intensity of 
the bumps behaves at infinity. There is, in fact, some competition among 
different properties of the potential, resulting in rich spectral possibilities. 

The prime numbers are straight related to both kind of potentials dis- 
cussed above, namely nonrandom and sparse. It is well known that the set 
of prime numbers is not finite and that there are arbitrarily large gaps be- 
tween consecutive primes, a property characterizing sparseness. By defining 
a potential Vo = 0,V^ = lifn+lisa prime number and V n = if not, 
and extending symmetrically to negative values of n, i.e., V_ n = V n , one gets 
an instance of finite valued sparse potential. The nonnegative values of the 
potential V are 01101010001010001 • ■ ■. 

The main aim of this Letter is to discuss numerically the spectral proper- 
ties of the corresponding discrete Schrodinger operator, which will be shortly 
referred to as "the prime Schrodinger operator." The rigorous analytical 
study of this system, although desirable, seems to be subtle and far from 
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trivial. Since in this Letter we consider this problem from just a numerical 
point of view, we find more appropriate to avoid the mathematical spectral 
terms, and instead shall try to use the physical ones localization (insulator) 
and derealization to characterize our results. 

By exploring dynamical as well as geometrical quantities related to the 
spectral properties of the model, we find a (rather smooth) transition from 
predominantly delocalized states to localized ones, as the potential intensity 
is increased. Nevertheless, some delocalized states always remain, no matter 
how large is the potential. Also, the combined use of geometrical and dy- 
namical tools, supportive and to some extent complementary to each other, 
seems to be applicable to other cases with coexisting spectral components. 

Here two different sets of numerical tools are used in such investigations. 
The dynamical tools are the average probability for the system to return to 
its initial state u = 6 n>no , concentrated on the site n , 

C{t) = \ f* \(u(s)\u )\ 2 ds, 



t Jo 

and the second moment 

d 2 (t)=J2(n-n ) 2 \u n (t)\\ 

n 

u = (u n ) ne -£. Recall that for large values of t it is expected that C{t) ~ t~ D ' 2 
and d 2 (t) ~ t 213 , with D 2 being the correlation dimension of the spectral 
measure associated with the initial state [14, 16 1. Continuous spectra are 



usually numerically characterized by < D 2 , /3; although it is possible to have 



exceptional systems with point spectrum and with any value < ft < 1 |L7 . 

Due to the rather irregular distribution of prime numbers |]J we suspect 
the quantum dynamics generated by Hy can strongly depend on the initial 
wave function, at least for small times; therefore, the initial condition was 
always concentrated on n = 0. These dynamical tools were restricted to 
0.5 < A, since for small A it is not possible to neglect logarithmic corrections 
to the algebraic behaviours of C(t) and d 2 {t) |J. 

The other kind of tools, geometrical tools, are the direct inspection of 
eigenf unctions and the Lyapunov exponent (LE) 7, also called inverse local- 
ization length, calculated for each eigenfunction of the finite basis approxima- 
tion. The typical basis size used here is 10 4 ; the robustness of the numerical 
calculations with respect to the basis size was verified in each case. Local- 
ization, i.e., point spectrum, is characterized by nonvanishing 7. 



From the Thouless formula H, 18 , in particular from its proof, one sees 
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that the LE at energy E can be estimated as 

1 N 

iV 3=1 

with Ej running over all eigenvalues (distinct from E) of the finite basis 
approximation of size N. 

In Fig. 1 it is shown the Lyapunov exponent 7 as function of the energy 
for some values of the potential intensity A. For A = 0.5 the vanishing of the 
LE for most energies E indicates the predominancy of delocalized states; by 
increasing the potential intensity the LE becomes different from zero, except 
in a left neighborhood of E — 0. This phenomenon was verified for all values 
of potential intensity considered (up to A = 40.0), i.e., the LE vanishes just 
below E — 0. This strongly indicates a transition from derealization to 
localization with a remaining set of delocalized states near E — 0. 

The above conclusions were also supported by numerical solutions of the 
Schrodinger equation and computation of dynamical quantities (also with 
lattice size up to 10 4 ). In Fig. 2 it is shown d 2 (t) for two values of A; d 2 (t) 
saturates for A = 6.0 (as well as C(t); not shown), indicating the presence of 
localized states; while for A = 0.5 the values of d 2 (t) increase until a critical 
probability density value is reached at the basis border (\u\ 2 < 10~ 6 was 
used), indicating the presence of delocalized states. Notice that, due to the 
presence of both localized and delocalized states for the prime Schrodinger 
operator, it can be difficult to retrieve information on its spectral type only 
from dynamical quantities. This justifies the choices of potential intensity 
values in Fig. 2, used to exemplify the spectral transition: for A = 0.5 most 
states are delocalized, while for A = 6.0 the energy band corresponding to 
delocalized states is very narrow. 

Of course a numerical vanishing LE is not exactly zero. The LE 7 p for a 
periodic case was computed with the above method; since it should be zero 
in the spectrum of the corresponding Schrodinger operator, any 7 < 7 P was 
in practice considered to indicate the presence of delocalized states for the 
prime case. This was used to estimate the length of the interval [— 6(A), 0] of 
the "delocalized band" as function of A. As another check for such transition, 
the dynamical quantities d 2 (t) and C(t) were computed for short times taking 
into account two different sets of eigenvectors: those whose eigenvalues were 
in the delocalized band [— 6(A), 0], and those in intervals with positive LE. 
The distinctions in their qualitative behaviours are evident. For example, in 
Fig. 3 the value A = 1.0 was fixed and d 2 was computed by restricting the 
dynamics to eigenvalues in the range [—0.6,0.0] (a subset of [— 6(1), 0]) and 
also to eigenvalues in the range [1.6,2.0] (corresponding to positive LE). In 
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the former case d,2(t) grows with time, while in the latter case it is clearly 
bounded; such results also support the presence of mobility edges in the 
spectrum of the prime Schrodinger operator. 

Fig. 4 shows some values of 6(A); the best fitting line is also shown and 
its slope is equal to —1. Therefore, it is numerically found that 6(A) ~ A -1 , 
at least for large A, so that it does not vanish and the spectrum of the prime 
Schrodinger operator should have a set of delocalized states for any poten- 
tial intensity. Notice that the eigenfunctions corresponding to eigenvalues in 
the range [— 6(A), 0] are extended over the finite bases, while they are expo- 
nentially localized if the corresponding LE is greater than zero (not shown 
here) . 

In summary, a delocalization-localization transition, with mobility edges, 
was found numerically for the one-dimensional prime Schrodinger operator, 
a transition as the intensity of the potential is increased; however, there re- 
mains a band of delocalized states whose length scales as ~ A" 1 , assuring 
the presence of delocalized states for any value of A. These results, and the 
calculations presented, suggest that the same approach could be applied to 
cases with other aperiodic potentials, for which one suspects to have simul- 
taneously continuous and point spectral components. 
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Figure Captions 



Figure 1. Lyapunov exponent 7 versus energy for the prime Schrodinger 
operator, a) A = 0.5; b) A = 1.0; c) A = 3.0; d) A = 6.0. 



Figure 2. The second moment d 2 as function of time for a) A = 0.5 and 
b) A = 6.0. Notice the different scales (for A = 6.0 the time evolution did 
not reach the basis border). 



Figure 3. The second moment d 2 as function of time for A = 1.0 (base 10 
log — log scale). The upper curve was obtained taking into account only the 
eigenvectors whose eigenvalues were in the range [—0.6,0.0]. Similar for the 
lower curve, but with eigenvalues in the range [1.6,2.0]. See also Fig. lb. 



Figure 4. Delocalized band length b as function of the potential intensity A 
(base 10 log — log scale). The best fitting line is also shown. 
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Figure 2 : a) , b) (top to bottom) 
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